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ABSTRACT 

The eigenfrequencies of nonradial oscillations are a powerful probe of a star's interior structure. This is 
especially true when there exist discontinuities such as at the neutron star (NS) ocean/crust boundary, as first 
noted by McDermott, Van Horn & Hansen. The interface mode associated with this boundary has subsequently 
been neglected in studies of stellar nonradial oscillations. We revisit this mode, investigating its properties both 
analytically and numerically for a simple NS envelope model. We find that it acts like a shallow surface ocean 
wave, but with a large radial displacement at the ocean/crust boundary due to flexing of the crust with shear 
modulus /j<P, the pressure. This displacement lowers the mode's frequency by a factor of ~ (/i/P) 1 / 2 ~ 0.1 
in comparison to a shallow surface wave frequency on a hard surface. The interface mode may be excited 
on accreting or bursting NSs and future work on nonradial oscillations should consider this mode. Our work 
also implies an additional mode on massive and/or cold white dwarfs with crystalline cores, which may have a 
frequency between the f-mode and g-modes, an otherwise empty part of the frequency domain. 
Subject headings: stars: neutron — stars: oscillations — white dwarfs 



1. INTRODUCTION 

A wide variety of nonradial oscillations have been previ- 
ously investigated for neutron stars (NSs) with the major- 
ity of these studies focusing on surface g-modes (oscillations 
whose restoring force is gravity; McDermott, Van Horn & 
Scholl 1983; Finn 1987; McDermott & Taam 1987; McDer- 
mott, Van Horn & Hansen 1988, hereafter MVH88; McDer- 
mott 1990; Strohmayer 1993; Bildsten & Cutler 1995, here- 
after BC95; Bildsten, Ushomirsky & Cutler 1996; Strohmayer 
& Lee 1996; Bildsten & Cumming 1998; Piro & Bildsten 
2004) because their frequencies are similar to many of the os- 
cillations observed from NSs (for example burst oscillations; 
Muno et al. 2001). These modes are confined to outer regions 
of the NS surface due to the shear modulus of the NS crust (as 
illustrated by BC95 and further discussed here), so that their 
frequencies are independent of the deep NS crust composition 
and small amounts of energy will lead to large amplitudes. 
However one mode has consistently been neglected, which is 
the wave analogous to a shallow water wave in the outermost 
parts of the star. Unlike a global surface wave (i.e. f-mode), 
we expect this wave to be confined to the NS ocean just like 
the g-modes, but a correct prediction of its frequency requires 
proper treatment of its eigenfunction at the ocean/crust bound- 
ary. 

A study of nonradial oscillations on non-rotating, non- 
magnetic NSs was presented by MVH88 for three-component 
NS models (consisting of an ocean, crust, and core), exam- 
ining g-modes, acoustic oscillations (p-modes), and f-modes. 
In addition, MVH88 reported a new set of oscillations, which 
they called "interface modes." These modes are unique in 
that they "live" (have amplitudes concentrated) around dis- 
continuities in the NS structure. For their three-component 
model, they found these modes at both the ocean/crust inter- 
face and the crust/core interface. We now undertake a de- 
tailed investigation of the ocean/crust interface mode so that 



it may be used for accreting NSs and in the other astrophys- 
ical context where a shear modulus is important, crystalline 
white dwarfs (WDs) (e.g. Hansen & Van Horn 1979, here- 
after HV79; Montgomery & Winget 1999). 

To focus our analysis we consider the interface mode for a 
two-component NS model consisting of an ocean and crust. 
We use this model because of the simplification provided by 
the plane parallel geometry of the outer NS layers. The most 
important property of the interface mode is how it couples 
with the top of the NS crust (Figure 1). MVH88 showed 
that the interface mode has a large and approximately con- 
stant transverse velocity in the ocean, similar to a shallow 
water wave. Indeed this mode can be thought of as a sur- 
face wave riding on a compressible crust. It is therefore not a 
toroidal mode as one might expect if the interface wave were 
like a shear wave (see Appendix A). As its pressure pertur- 
bation rolls over the crust's surface, the crust is compressed 
and sinks. This flexing causes the interface mode to have a 
large, negative radial displacement at the discontinuity (see 
the right diagram of Figure 1) in contrast with the typical 
boundary condition for surface water waves, which treats the 
crustal boundary as solid with no radial displacement (see the 
left diagram of Figure 1). Since the pressure perturbation also 
corresponds to a small, positive swell at the top of the ocean, 
the radial displacement eigenfunction must necessarily have a 
single node within the ocean. 

This simple picture describes the interface mode quite well. 
The frequency depends on the properties at the base of the 
ocean, like a surface ocean wave, but with a frequency that is 
altered by the non-zero radial displacement at the ocean/crust 
boundary. This boundary condition reduces the frequency by 
~ (jjl/P) 1 ' 2 ~ 0.1 from what is expected if the crust were 
"solid," where fi is the shear modulus and P is the pressure 
at the top of the crust. Most of the energy of the mode is in 
the ocean, just as in a shallow surface wave. 




FIG. 1 . — A schematic diagram of an ocean surface wave showing how it is altered by the flexing of the crust. The vectors show the wave velocity. In both 
cases the transverse velocity is approximately constant within the ocean and dies off quickly into the crust. If the crust is treated as if it is "solid" the radial 
displacement goes to zero at the ocean floor (left diagram). As shown on the right, the interface mode has a large, negative radial displacement at the crust which 
is actually larger than the corresponding swell near the ocean surface. This implies that its radial eigenfunction must necessarily have a single node within the 
ocean. 



We derive the adiabatic perturbation equations for a plane 
parallel layer in hydrostatic balance with a shear stress in §2. 
These are used in §3 to analytically estimate the interface 
mode's radial eigenfunction and frequency. We follow this 
with a numerical study in §4 that confirms our analytic results 
and also provides more detail about the eigenfunctions. We 
conclude in §5 with a discussion of possible uses for the in- 
terface mode in astrophysical contexts. Maybe most exciting 
among these is an additional mode of oscillation in WDs with 
crystalline cores. 

2. ADIABATIC PERTURBATION EQUATIONS 

Since the pressure scale height in the NS ocean, h — 
P/{ps) ~ (2—10) x 10 2 cm, is always much less the NS ra- 
dius, R w 10 6 cm, we approximate the surface as having con- 
stant gravitational acceleration, g = GM/R 2 (neglecting gen- 
eral relativity), and plane parallel geometry. We use z as our 
radial coordinate andx as the transverse coordinate. An addi- 
tional useful coordinate in this geometry is the column depth, 
denoted as y (defined by dy — —pdz and with cgs units of 
g cm" 2 ). 

The interface mode's period is set by physics at the 
ocean/crust boundary, and is insensitive to the outer bound- 
ary condition. At these depths the thermal timescale is longer 
than the mode period (since ?th ~ c p yT /F and y is large; where 
c p is the specific heat capacity and F is the flux) so that the 
perturbations are adiabatic. The governing equations of the 
envelope are conservation of mass, dp/dt + V ■ (pv) = 0, and 
momentum, 
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where ct, 7 is the stress tensor, and z is a radial unit vector (re- 
peated indices imply the Einstein summation convention un- 
less noted otherwise). We neglect the NS spin so that we can 
focus on the effects of the crust. Rotation can be simply in- 
cluded using the "traditional approximation" (e.g. Bildsten, 
Ushomirsky & Cutler 1996). 



We make Eulerian perturbations of the conservation equa- 
tions, substituting for the dependent variables Q — ► go + SQ, 
where Qo denotes the static background (this subscript is 
dropped in subsequent expressions). Perturbations are as- 
sumed to have the form SQ = 6Q(z)cxp(ik x x + iuit), where 
k 2 — 1(1 + l)/R 2 is the transverse wavenumber and u> is the 
mode frequency. These Eulerian perturbations are related to 
Lagrangian perturbations by AQ = SQ + £ ■ VQ. Since the 
background model has no fluid motion we substitute v — > <5v = 
d£/dt = iu>£, where £ is the local Lagrangian displacement. 
We use the Cowling approximation and neglect perturbations 
of the gravitational potential (Cowling 1941), an excellent ap- 
proximation in the thin outer layers of the NS. Keeping terms 
of linear order, we find 

Sp + V-(p$)=0, (2) 

for continuity. 

In making linear perturbations of equation (1) we include 
the effects of a finite shear modulus (such as in the NS crust). 
Relative displacements then introduce a shear component to 
the Lagrangian stress tensor given by (Landau & Lifshitz 
1970). 



(3) 



where p is the crystalline shear modulus as calculated by 
Strohmayer et al. (1991) for a classical one-component 
plasma (OCP), 

0.1194 n,(Ze) 2 



P : 



(4) 



l+0.595(173/r) 2 a ' 

where m is the ion number density, a = (3/4nm) 1 / 3 is the 
average ion spacing, Z is the charge per ion, and 

~ ak B T T 8 /4 1.30 / \ A 1 1 1 



v 10 9 g cm~ 3 / 

is the dimensionless parameter that determines the liq- 
uid/solid transition for the OCP, where k B is Boltzmann's con- 
stant, Tg = r/10 8 K, and A is the number of nucleons per ion. 
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This transition occurs at T w 173 (Farouki & Hamaguchi 1993 
and references therein, following the work of Brush, Sahkin 
& Teller 1966). Since the pressure in the crust is dominated 
by degenerate electrons we rewrite p as 

-2 /z x2/3 



and 



p 



1.4 x 10" 



l+0.595(173/r) 2 V30, 

In the crust p/Pis fairly independent of temperature (except 
for a small dependence due to the factor of T in the denom- 
inator), so we typically substitute = fJ-/P an d assume that 
po is constant with depth. 

In the adiabatic limit, AP/P = Ap/(Y\p), where T\ = 
(d\nP/d\np) s is the adiabatic exponent. The Lagrangian 
stress tensor is 

Aa ij =PT l (V-Z)5 ij + Ao s ij , (7) 
and the Eulerian stress tensor is 

6(T U = Aai j -pg£ z 6 u , 
so that the perturbed momentum balance equation is 
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(9) 



Equations (2) and (9) describe the nonradial oscillations. 

The fluid viscosity in the ocean barely resists shear so that 
the Lagrangian stress tensor has only diagonal components, 
Aoij = -APSjj. Breaking equation (9) into radial and trans- 
verse components, 



■ pu\ - 
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(10) 



- P u?i x =-ik5P. (ii) 

Combining equations (2), (10), and (11) results in the stan- 
dard nonradial oscillation equations for an inviscid fluid in 
hydrostatic balance and plane parallel geometry, 



dz 



d 5P 
~dz~~P 



1 



Tih 
1 5P 

hT 



\ LU 2 



1 \ SP 



where 



N 2 = -g 



dlogp 1 dlogP 





P ' 






8 / 









(12) 



1 £ , (13) 



(14) 



dz Ti dz 

is the Brunt-Vaisala frequency which measures the internal 
buoyancy of the envelope. Normal modes of oscillation are 
found by assuming AP = at the top boundary (a fluid ele- 
ment riding at the top feels no pressure perturbation) and then 
shooting for the bottom boundary condition. This top con- 
dition, though not unique, is fairly robust since little mode 
energy resides in the low density upper altitudes of the ocean 
(BC95). 

When pj^O, expanding equation (9) in both the radial and 
transverse directions, and using the adiabatic condition and 
equation (2) to eliminate Sp and SP, results in 
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the crustal mode equations derived by BC95. 

The second derivatives in equations (15) and (16) imply 
additional boundary conditions at the ocean/crust interface. 
We require that £ z is continuous, so that there is no "space" 
between the two layers, and that Aa zz is continuous, other- 
wise there is an infinite radial acceleration across the bound- 
ary. Furthermore, Aa xz = because the ocean, with no shear 
modulus, cannot sustain finite transverse shear. We also re- 
quire the additional boundary conditions that £ z k(,k0 at 
some depth deep within the crust. Though this is not neces- 
sarily the case, it is helpful for estimating the interface mode's 
frequencies. 

3. ANALYTIC ESTIMATES 

We now use these differential equations to analytically de- 
rive the frequencies and eigenfunctions of the interface mode. 
As we explain in Appendix A, the interface wave is not well 
described by WKB analysis, so we proceed with use of analo- 
gies to shallow ocean surface waves. In §3.1, we solve equa- 
tions (12) and (13) for a surface wave, and consider how its 
frequency changes when £ z ^ at the ocean floor (as sug- 
gested by Figure 1). In §3.2 we solve for the radial crustal 
eigenfunction, which we use in §3.3 to connect to the ocean 
perturbation and find the radial displacement at the interface. 
In §3.4 we summarize our frequency prediction for the inter- 
face wave and show how it depends on the temperature and 
composition of the NS crust. 

3.1. Surface Wave Analytics 

We assume that the surface wave in the ocean has a constant 
transverse displacement so that £ x = £ xX and d£/dz = (where 
we use the subscript "t" to denote the "top" of the ocean). At 
the top the boundary condition is AP = = SP — i z pg, and 
by combining this with transverse momentum conservation, 
equation (1 1), we relate the surface radial and transverse dis- 
placements 



u> 



£z,t — Tj—£,x,t- 
ikg 



(17) 



Using equation (11) to eliminate SP from equation (12), and 
then eliminating £t )t using equation (17), 



dz 



£z,t 
h 



8hk 2 x 



1 



(18) 



The electrons in the deep NS ocean are degenerate and rela- 
tivistic so that the equation of state can be approximated by 
a P = Kp 4 / 3 polytrope, where K is a function of the mean 
molecular weight per electron. Using hydrostatic balance, the 
pressure scale height is 



ft(z) = i(z t -z), 



(19) 



1 /3 

where z t = 4-Kp c ' /g, the subscript "c" denotes the top of the 
"crust", and z = at the ocean/crust interface (we also make 
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the approximation that p c ^> p t ). Substituting this for h in 
equation (18) and integrating gives, 

where the constant of integration is set to assure that £ z = £ zt 
atz = z,. 

The frequency of this surface wave is set by specifying the 
boundary condition £ Z:C = £ z (z = 0). Substituting this into 
equation (20) we find a frequency 

tJ _ (4/r 1 + i)-'g*& (21) 

1 — £z,c/£z,t 

To understand this result it is helpful to compare it to the fre- 
quency for a mode with £ z , c = as expected for a solid floor, 
which we denote as ujq. This results in 

-l 



w n = 



■ + 1 



gk 2 x z t - 



(22) 



For =4/3 we find that = gh c k 2 , where ft c is the scale 
height at the bottom of the ocean, just as expected for the 
shallow surface wave. The interface mode's frequency is then 



og/(l-W&,t)- 



(23) 



We expect &,c/&,t < because a positive pressure perturba- 
tion will compress the crust (see the discussion in § 1 and Fig- 
ure 1). Equation (23) shows that flexing of the crust results in 
a lower frequency than what would be expected if the ocean 
floor were solid. 

3.2. The Radial Crustal Eigenfunction 

To simplify equations (15), so that it can be solved analyti- 
cally, we order the terms and drop those which are small. We 
first relate the ordering of the transverse and radial displace- 
ments within the crust. The ocean and crust slip with respect 
to one another so that the transverse displacement changes 
discontinuously at the boundary. We therefore set ^ c = A£ tt 
where A is a dimensionless constant which we call the "dis- 
continuity eigenvalue." It determines how £ x changes at the 
interface, and is an eigenvalue set by the bottom boundary 
condition within the crust. We then rewrite equation (17) as 



f - igkx \f 

Sx,c — — — A ?; I • 



(24) 



At the top of the crust, the displacements must satisfy Aa xz 
so that 



dz 



(25) 



where we have used /io = n/P- Assuming that 0(d£ x /dz) ~ 
0(£x C /h) within the crust, we combine equations (24) and 
(25) to find 



0(A) ~ O 



(26) 



The largest reasonable frequency is that of a discontinuity 
g-mode at the boundary (cu 2 w ghk 2 ) for which the terms in 
equation (15) with J 1 appearing are negligible. Using order 
of magnitude estimates of k x h ~ 10~ 4 - 10 -3 and p Q « 10~ 2 , 
from equation (6), we identify the two highest order terms of 
equation (15), resulting in 



dz 2 



h dz ' 



(27) 



as the simplified radial differential equation. 

Since Ti w 4/3, we introduce a new radial variable, s, de- 
fined by s = Ah which is set equal to s c at the ocean/crust 
interface and increases going down so that ds = —dz. Integra- 
tion of equation (27) results in two constants of integration. 
The first is set by requiring that £ z = £ ZjC at s = s c and the sec- 
ond by requiring that £ z = at some bottom depth, which we 
denote as s\» so that 



(s b /s)' 
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(Sb/Sc 



1 



(28) 



If we place the lower boundary condition infinitely deep 
within the crust, so that s\, >• s c , equation (28) simplifies to 



(29) 



identical to what BC95 estimate for the crustal eigenfunction 
(in their case in the context of a g-mode extending into the 
crust). 

3.3. Connecting the Modes at the Ocean/Crust Boundary 

We solve for £ ZjC , and estimate the mode frequency using 
equation (21), by demanding that Act,, be continuous across 
the ocean/crust boundary to avoid infinite radial accelerations. 
Expanding this component of the stress tensor, 

^=(r,-^)(V.O + 2Mof. (30) 

We evaluate this condition using the eigenfunctions we have 
found for both the ocean and crust. In the ocean, above the 
interface, 



Act, 



above 



— T~ (£z,c — £z,t)- 

h c 



(31) 



In the crust we use equation (28) to evaluate d^Jdz. The 
term involving ik x £_ x can be ignored because it is order (k x h c ) 2 
less than the d£ z /dz term (see Appendix B), so that below the 
interface, 



Act,. 



below 



3r 1 /4 + /x 

1 - (s c /s b y 



h c 



(32) 



where we have substituted s c = 4h c . Setting equations (31) 
and (32) equal, we solve for the radial displacement at the 
crust, 



3^/4 + ^ 



1 



1 

Mo' 



(33) 



.l-(*cAb) 3 

where the last estimate is made for Ti =4/3 and s\, » s c . 

3.4. Interface Wave Frequency Estimate 

From our estimate for the eigenfrequency, equation (23), 
along with £ z , c /£ z ,t ~ — I/Mo. equation (33), we find the sim- 
ple result that 



to 



2 



(34) 



i li 

The frequency of the interface mode is a factor of /U ' ^0.1 
less than a shallow surface wave in the NS ocean. 

It is an interesting exercise to use this result to predict the 
interface wave's frequency and dependence on the ocean/crust 
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boundary's properties. At this depth, the pressure is dom- 
inated by relativistic and degenerate electrons with a Fermi 
energy of (Ep > m e c 2 » k B T) 



E F = 4.l MeV 

Since the pressure is P ■ 
is 

h c = 2620 cm 



2Z 
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(35) 



10 9 g cm- 
n e Ep/4, the scale height at the crust 
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(36) 



Using ujq = gh c k 2 and substituting Y from equation (5) to 
eliminate p, we find the interface wave has a frequency of 
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(37) 



For an accreting NS with a relatively cold ocean (similar to 
BC95 or what we consider in §4.1), this frequency is between 
the g-mode frequencies, < 6 Hz (BC95), and the f-mode fre- 
quency, ~ 10 3 Hz, an otherwise empty part of the frequency 
domain. Furthermore, the interface mode frequency is inde- 
pendent of Z and has a simple scaling with both A and T, 
making this mode a useful diagnostic for measuring the prop- 
erties of the NS crust. 

4. NUMERICAL SOLUTIONS 

We now compute eigenfunctions and frequencies for the in- 
terface mode using numerical integrations of the mode equa- 
tions on a two-component, ocean/crust NS model. In §4.2 we 
show that these compare very well with the analytic work of 
§3. 

4. 1 . NS Ocean and Crust 

We assume that the NS has a mass of M = 1 .4M and radius 
of R = 10 km. The two components which make up our NS 
envelope are a liquid ocean and a degenerate crystalline crust 
(see Figure 1). For accreting NSs, the composition of each 
of these components is determined by nuclear burning pro- 
cesses which occur on the surface, especially the rp-process 
during type I X-rays bursts (Wallace & Woosely 1981; Schatz 
et al. 1999, 2002) and the hotter burning during superbursts 
(Cumming & Bildsten 2001; Schatz, Bildsten & Cumming 
2003). Woosley et al. (2004) follow the accretion and sub- 
sequent unstable burning over multiple type I bursts in nu- 
merical simulations. Estimating their results for an accreted 
composition of 0.05 solar metallicity and an accretion rate of 
1.75 x 10~ 9 M Q yr" 1 (which they denote model "zM") after 
14 bursts, we assume a composition with mass fractions of 
64 Zn(0.56), 68 Ge(0.35), 104 Pd(0.07) and 106 Cd(0.02), for both 
the ocean and crust. At sufficient pressure and density, this 
material will change from liquid to solid. This transition oc- 
curs at T w 173 (see §2.1) which implies a critical density of 

Pc =2.5x 10 9 gcm- 3 (j 



64 



r 

173 



(38) 



for an OCR We assume that this is a good enough approxi- 
mation even though we use a multi-component composition. 



Due to the strong dependence of this critical density on Z, 
in cases where burning results in ashes of especially heavy 
elements (such as the rp-process at high accretion rates or 
during superbursts) crystallization can happen at significantly 
shallower depths. This decreases the energy needed for large 
mode amplitudes. 

The envelope profile is described by the equation of radia- 
tive transfer, which after substituting y is 



F = 



4acT 3 dT 
3k dy ' 



(39) 



where a is the radiation constant and k is the opacity. The 
opacity is set using electron-scattering (Paczyhski 1983), free- 
free (Clayton 1993 with the Gaunt factor of Schatz et al. 
1999), and conductive opacities (Schatz et al. 1999 using the 
basic form of Yakovlev & Urpin 1980) in the liquid region. In 
the crust we make no distinction and use the same conductive 
opacity. 

Equation (39) is integrated starting at an outer column depth 
(y w 10 7 g cm~ 2 ) and into the envelope. The outer tem- 
perature is set near 10 8 K, a value for which the solutions 
are minimally sensitive because of their radiative zero nature 
(Schwarzschild 1958, or see the discussion in Piro & Bild- 
sten 2004). Previous studies of constantly accreting NSs have 
shown that the interior thermal balance is set by electron cap- 
tures, neutron emissions, and pycnonuclear reaction in the in- 
ner crust (Miralda-Escude, Paczyhski & Haensel 1990; Bild- 
sten & Brown 1997; Brown & Bildsten 1998) which release 
« 1 MeV/m p « 10 18 erg g' 1 (Haensel & Zudunik 1990). De- 
pending on the accretion rate and thermal structure of the 
crust, this energy will either be conducted into the core or 
released into the ocean such that for an Eddington accretion 
rate up to w 92% of the energy is lost to the core and exists 
as neutrinos (Brown 2000). We therefore set the flux in the 
ocean to a fiducial value of 10 21 erg s~'cm~ 2 (as expected 
for an accretion rate of <~ 1O~ 9 M yr _1 , about one-tenth the 
Eddington rate). We solve for p using the analytic equation 
of state from Paczyhski (1983) which is applicable to both 
the ocean and crust. The profiles are continuous across the 
ocean/crust interface because we assume the same conductive 
opacity in each region. Figure 2 shows the envelope profile for 
the variables most important for determining the mode prop- 
erties. The dashed line denotes the location of the crust (which 
we assume is at a column depth of y c 
density of p c w 10 9 g cm~ 3 ). 



10 13 g cm 2 , near a 



-2 



4.2. Numerical Integrations and Comparisons with 
Analytics 

Using the envelope profile of §4.1, we numerically solve 
for the interface wave eigenfunctions and eigenfrequencies. 
We treat each integration as a two eigenvalue shooting prob- 
lem, first choosing u and A and then shooting for the bottom 
boundary condition. We integrate equations (12) and (13), be- 
ginning at the top of the ocean at a column of y t « 10 8 g cm -2 
(where the unstable burning of X-ray bursts typically occurs; 
Bildsten 1998) where we set AP — 0. The boundary condi- 
tions at the ocean/crust interface are set by requiring £ z to be 
continuous, Aa xz = (which sets d^ x /dz), and Aa zz continu- 
ous across the boundary (which sets d£ z /dz). The discontinu- 
ity of £x is set using A. We then integrate equations (15) and 
(16) into the NS crust (we assume that dpo/dz ~ dT\/dz ~ 
and simply set po = 10~ 2 for these examples) down to a depth 
within the crust where we require both £ A and £ z to be zero. We 
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FIG. 2. — Temperature, density, and Brunt-Vaisala frequency for our two- 
component neutron star envelope. The dashed line denotes the depth of the 
ocean/crust boundary, which lacks a discontinuity because we use the same 
conductive opacity in each region. The flux is assumed to be constant at 
F = 10 21 erg s~ 1 cm" 2 with no other heat sources. This is a fiducial value 
expected for accreting neutron stars which have their internal thermal bal- 
ance set by electron captures, neutron emissions, and pycnonuclear reac- 
tions in the inner crust. These reactions release a 1 MeV/m p with typi- 
cally f» 0. 1 MeV / nip making its way into the upper crust (Brown 2000). The 
power law of the density (p oc }> 3 / 4 ) shows that the pressure in the crust is set 
by degenerate, relativistic electrons. 



leave the bottom boundary depth as a free parameter and study 
how the eigenvalues and eigenfunctions change with this vari- 
able. 

In Figure 3 we show example eigenfunctions with a bot- 
tom boundary at y\, — 10 16 g cm" 2 . On the same plot we 
also compare the analytic eigenfunctions from §3.2 and Ap- 
pendix B. These use the numeric envelope model to set the 
local scale height and T\ = 1.338 (its value at the ocean/crust 
boundary) to set the eigenvalues. The numerical eigenval- 
ues are w/(2tt) = 23.43. Hz and A = -1.5 x 1(T 4 , which 
are close to the analytic values of oj/(2tt) — 24.46 Hz and 
A = —1.8 x 10~ 4 . As shown in Figure 1, the radial displace- 
ment is found to be very large at the interface, relative to its 
value at the ocean surface. The energy density (per logarithm 
of pressure) is 



(40) 



where £ 2 = |£ x | 2 + |£.| 2 , which shows where the mode "lives." 
The energy density is also strongly concentrated at the inter- 
face, consistent with the result that the mode's properties are 
set by the envelope at this depth. 

As a comparison, we plot £ z for both the correct and the 
"solid" ocean/crustal boundary conditions in Figure 4. For 
the latter we find a frequency of uq/(2-k) = 179.62 Hz. The 
linear plot of the amplitudes shows how extreme the radial 
displacement is at the crust for the interface mode. We omit 
comparing £ x and the energy density for the modes since they 
are almost identical in the ocean. 




^0 ^r- 

o 35 



Column Density (g ci 



FIG. 3. — A comparison of the numerical eigenfunctions (solid curves) 
and the analytic approximations (dashed lines) for the interface mode us- 
ing the envelope shown in Figure 2, with the bottom boundary set at 
yt, = 10 16 g cm" 2 . The numerical solutions find uj/(2tt) = 23.43 Hz and 
A = —1.5 X 10" 4 , while the analytic values are u)/(2ir) = 24.46 Hz and 
A = —1.8 X 10" 4 , using equations (21) and (B7), respectively. The ana- 
lytic solutions use the numeric envelope to set the local scale height, and 
T[ = 1.369 to set the eigenvalues (this is the value of Ti at the base of the 
ocean). The radial eigenfunction, £ z , shows a single node in the ocean and a 
large displacement at the crust (y c = 10 13 g cm" 2 ), consistent with the dis- 
cussion in § 1 . The energy density (bottom panel) shows that the mode energy 
is concentrated in the ocean. 
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FIG. 4. — A comparison of the £ z eigenfunction between a correctly calcu- 
lated interface wave (solid curve) and the case of when the crustal boundary 
is set to £ z = (a standard shallow ocean wave; dashed curve). The latter case 
has an eigenfrequency of o;o/(27r) = 179.62 Hz and no node in the ocean, in 
contrast to the interface wave's eigenfrequency of ui/(2ir) = 23.43 Hz. The 
linear plot of the displacements (bottom panel) emphasizes the large displace- 
ment at the boundary for the interface wave (note the other mode is magnified 
by a factor of 20). 
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FIG. 5. — The properties of the interface mode as a function of the 
bottom boundary depth, using the envelope of Figure 2 with the crust at 
y c = 10 13 g cm -2 . The solid curves show the analytic estimates, using equa- 
tions (21), (B7) and (33) (from top to bottom), while the points mark the 
numerical calculations. On the far left we find what is expected for a shal- 
low ocean wave, and as the curves move right (to larger depths) they become 
more like an interface wave. 



To study the dependence of the mode's properties on the 
bottom boundary condition within the crust, we compare 
the numeric and analytic results as a function of the bottom 
boundary depth, y\» in Figure 5. On the left hand side is the 
limit when there is no penetration into the crust, in this case 
both the ratio £ z , c /£z.t and A are zero and uj/(2tt) ~ 200 Hz, 
the values which would be approximated for a shallow surface 
wave (see Figure 4). As the bottom boundary is set deeper, the 
frequency decreases, becoming nearly, but not quite equal to 

- [ij 2 x 200 Hz ~ 20 Hz. 

5. DISCUSSION AND CONCLUSIONS 

We have discussed the ocean/crust interface wave, inves- 
tigating its properties both analytically and numerically for 
a two-component NS envelope. We show that it can be un- 
derstood as a shallow surface wave, with a different radial 
displacement boundary condition at the bottom of the ocean. 

1 12 

This reduces the frequency by roughly <~ fj,J ^0.1 from 
what would be expected if the crust were instead "solid." In 
general the interface wave is a global mode of oscillation, so 
that the entire star must be taken into account to calculate the 
mode to high precision (see discussion in Appendix C). For- 
tunately, by using our analysis, a fairly good estimate of the 
mode frequency can be made just from considering the NS 
properties right near the ocean/crust interface because the ma- 
jority of mode energy (> 99%) resides in the NS ocean. 

The interface wave should be included in future studies of 
nonradial oscillations on both accreting and bursting NSs. 
For a non-rotating NS this would result in frequencies of 
~ 10 — 20 Hz which would then be altered and split due to 



the NS spin (Bildsten, Ushomirsky & Cutler 1996). The ob- 
servation of multiple mode periods would greatly assist our 
understanding of the NS crust. The g-modes are sensitive to 
the depth of the crust (since it provides a bottom boundary), 
while the interface wave is sensitive to both the depth and the 
value of the shear modulus. This provides information about 
the temperature and composition of the crust, which is de- 
termined by unstable nuclear burning processes (X-ray bursts 
and superbursts) that happen in the ocean. 

Since the example considered in this work is a relatively 
"cool" envelope, the g-mode frequencies (< 6 Hz) are well 
separated from the interface wave frequency (~ 10 — 20 Hz). 
On a bursting NS, which has temperatures in the upper at- 
mosphere as high as 10 9 K, at least one g-mode has a fre- 
quency higher than the interface wave. The interface wave 
would then have an additional node in its eigenfunction, due 
to the relation between g-mode frequency ordering and num- 
ber of nodes, but its frequency would be unchanged as long 
as the ocean stays at a temperature of ~ 4 x 10 8 K (since this 
is where the mode's energy is concentrated). Woosley et al. 
(2004) show that the thermal wave during an X-ray burst does 
not penetrate deeply into the ocean, so that our estimates for 
the interface mode frequency still hold. 

The other important astrophysical context where our work 
can be applied is for old and/or massive WDs that have crys- 
talline cores. This may be seen in the pulsations of BPM 
37093 (Kanaan et al. 1992), a high-mass ZZ Ceti star that 
should be considerably crystallized (Winget et al. 1997). 
There is also the exciting possibility that some of the pulsators 
discovered through the Sloan Digital Sky Survey (Mukadam 
et al. 2004) may be massive enough to be crystallized. Scaling 
equation (37) to the fiducial interior properties of a w 1 . 1M Q 
WD, we find a crustal mode period of 



Pw95 s 



173 

^ 1/2 
16 



1/2 



10 7 K 
T 

R 



1/2 



5x 10 8 cm/ [1(1 + 1) 



1/2 



(41) 



which for the I = 1 mode is between the periods of the f-mode 
(~ 10 - 20 s) and g-modes (> 300 s). HV79 (and more re- 
cently Montgomery & Winget 1999) previously calculated the 
spectrum of g-modes and torsional modes expected for these 
stars, concentrating on the effects of crystallization. HV79 
do not single out any one mode as a "crustal mode," but in 
their crystallized models their lowest order g-mode (which 
they denote -iSf) matches many of the properties we iden- 
tify as being characteristic of a crustal mode. This includes 
a large radial displacement at the crust and a nearly constant 
transverse displacement in the WD "ocean" (see Figures 1 and 
2 of HV79). Furthermore, HV79 find the period of their _ t 5/ 
mode to be between the f-mode and the other g-mode periods, 
in confirmation of our above discussion. 

We thank both Phil Arras and Omer Blaes for helpful dis- 
cussions. We also thank the referee, Michael Montgomery, 
for pointing out the HV79 paper and its relevance to our work. 
This work was supported by the National Science Foundation 
under grants PHY99-07949 and AST02-05956, and by the 
Joint Institute for Nuclear Astrophysics through NSF grant 
PHY02-16783. 
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APPENDIX 



A. INFLUENCE OF A FINITE SHEAR MODULUS 

The interface wave's properties are set by how it is influenced by the finite shear modulus of the NS crust. Unlike g-modes, 
which see the crust as solid, the interface wave reacts to the crust as if it is a flexible and compressible surface. We now explain 
why each of these cases occur. 

Using /io = [i/P, and assuming that dfio/dz « dT\/dz « 0, equations (15) and (16) can be rewritten as 



dz 2 



3M0 



h 



+ik x £ x 



(1-Ti) 2/zp- 
h 3 h 



-ik x 



dz 



(f +ri ) 



(Al) 



and 



-Mo 



dz 1 



, < 1 Mo 1 



Tik 2 , 



Mo^u 
h dz 



- ik 



di z /Mo . r 

& lT + ri 



(A2) 



The simplest analysis of these equations arises when perturbations scale like ~ exp(/A: z z) and dispersion relations are in the WKB 
limit (when the radial wavenumber k z ^> l/h). This results in toroidal and spheroidal modes with frequencies 



^shear " 



w s 2 ound = (ri+4Mo/3)(P/p)fe 2 



(A3) 
(A4) 



where k 2 = k 2 +k 2 is the total wavenumber. The first is a shear mode (with a frequency of ~ 10 2 — 10 3 Hz) and the second is a 
sound wave slightly modified by 4/i /3 (~ 10 4 - 10 5 Hz). These imply that a g-mode with a frequency of < 6 Hz (BC95) will 
be reflected at the interface because of the frequency mismatch. Unfortunately, this analysis cannot be extended to the interface 
mode because it has k z < 1 /h, in contradiction to the WKB limit. 

To see that a mode with a constant transverse velocity (like the interface wave) is not excluded from the crust we compare the 
transverse momentum terms from equation (9). The shear modulus term has two pieces, the larger one is d6a xz /dz ~ [id £, x /dz 2 . 
This term strongly affects the mode's eigenfunctions when it is comparable to the transverse acceleration of the mode, pu 2 t; x . 
Using d/dz~k z this implies a critical shear modulus of n cr it/P ^oj 2 / (ghk 2 ). In the case of g-modes k z ~ l/h andw 2 K.AN 2 {hk x ) 2 
so that the critical modulus is (BC95) 



Merit 



(g-mode) 



6l(l + l)k B T (h 
Z \R 



Substituting values near the ocean/crust interface for degenerate, relativistic electrons, 



Merit 



(g-mode) - 10-'% 2 1(1 + 1) 



64 



30\ 



2/3 



r 

173 



(A5) 



(A6) 



which is clearly less than the value of [ijP from equation (6), implying that the g-modes will be excluded from the crust. 

One the other hand, the interface wave has a constant transverse velocity so that k z ~ \/R. Using J 2 w [iQghk 2 , the interface 
wave has a critical shear modulus of 



Merit 



(interface wave) ~ + 1), 



(A7) 



which is just the actual shear modulus present at the boundary, so the wave must penetrate the crust. This shows that in the 
interface wave the transverse momentum is perfectly balanced by the shear stress, so that the ocean/crust boundary behaves like 
a flexible and compressible surface. 



B. TRANSVERSE CRUSTAL EIGENFUNCTION AND DISCONTINUITY EIGENVALUE 

To estimate the transverse crustal eigenfunction, we order the terms from equation (A2) and drop those which are small. Using 
equations (24) and (26) we first relate the ordering of the transverse and radial displacements within the curst. Dividing equation 
(A2) through by £ z , multiplying through by h 2 and using the above relations, the ordering of each of the terms are 

o( PQ k x h) ~ o{uj 2 h 2 k x /g) + o(M> 3 ) + o(k 3 x h 3 ) + o{k x h) 

+O( f i k x h) + 0(k x h) + 0(^ Q k x h) + 0{k x h). (Bl) 
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The size of these quantities are approximated using k x h ~ 10~ 4 — 10~ 3 and /j,q w 10~ 2 . Just as with the radial equation considered 
in §3.2, the terms which contain ui 2 are negligible (as is seen from considering u 2 anywhere in the range of ~ fioghk 2 — ghk 2 ). 
To find a relation which includes ^ x we must consider both first- and second-order terms of equation (A2) to get 



dz 2 h dz Z \h h 



(B2) 



as a new simplified differential equation. 

We substitute the eigenfunction for £ z , equation (28) into equation (B2) to solve for £ x . The two constants of integration are set 
by imposing that Aa xz = at the top of the crust and that £ x (just like £ z ) goes to zero at s^,, so that 



(sb/sc) 3 



- 1 



+ 



+7 



7 1 1 



3/xo Sb 3/io V^b 
3/i + 4-3ri 



where we have introduced dimensionless constants that only depend on the star, 

2 / u + 4-3r 1 .y c 1 
a = ■ 



/? = - 
7=1 



2/<o 



1 

Mo' 

Substituting s = s c and dividing by equation (17) we find 

^Wg gz.c f 

" (v*c) 3 - 1 e, t I 



(B3) 

(B4) 
(B5) 
(B6) 



* I -1 



■0 



* I 1 



+7 



1_* 



(B7) 



for the discontinuity eigenvalue. 

We use equation (B3) (or consider the infinite crust limit of Appendix C), to see that the size of £t jC /£z,c ~ k x s c . This shows 
that when we consider terms in equation (30), the £ x terms are of order 



(*A) 2 «i, 



(B8) 



d£ z /dz 

smaller than the £ z terms and can thus be ignored. 

C. INFINITE CRUST LIMIT 

Just as we did for the radial eigenfunction for equation (29), we take the s\, ^> s c limit on equation (B3) to find the transverse 
eigenfunctions in the infinite crust limit, 



' ikx^C^Z: 



2/zo + 4-3ri /j c \3 



3/xo ^ s 
3jxo_+4-3Ti /j c \2' 
2/io V s J 

3 



(*)'-§(# 



(CI) 



which corresponds to 



A = — 



cj 2 s c ^, c / 5Mo+4-3r 1 

8 &,t V 6/xo 
5 u) 2 s c 
'6 g 



(C2) 
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for the discontinuity eigenvalue. 

We use these eigenfunctions to estimate the mode's energy density, dE/dlogy, from equation (37). We consider the limit of 
j^ooto see if the energy density falls off with increasing depth. If this happens, it implies a natural depth at which to set a 
lower boundary and provides confidence that the interface mode is determined by conditions local to the ocean/crust boundary. 
Comparing equations (29) and (CI) we find that at depths of order the NS radius £ x will dominate so that £ 2 oc 1 js 4 '. For fully 
degenerate and relativistic material P oc y oc p 4 / 3 . Using h = P/{pg) and h oc s we find that y oc s 4 . This means that at the 
greatest possible depths dE/d\ogy ~ £ 2 y ~ constant, so that, unfortunately, deep regions of the NS much be taken into account 
to calculate the interface mode to very high accuracy. In practice, this is not necessary because there are structural changes within 
a realistic NS crust that introduce new bottom boundaries. Due to the large amount of energy concentrated in the ocean (> 99%, 
see Figure 3), our frequency estimates are still good even with this difficulty. 
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